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Abstract 

We present the results of a quenched lattice QCD study of the coupling g, in the 
static heavy quark limit. After combining this with our previous results obtained by 
using propagating heavy quarks with a mass around the physical charm quark, we 
are able to interpolate to the 6-quark sector. Our results are = 0.48 ±0.03 ±0.11, 
g b = 0.58 ± 0.06 ± 0.10, and g B . Bw = 47 ± 5 ± 8. 



PACS: 12.38.Gc (Lattice QCD calculations), 13.75. Lb (Meson-meson interactions) 



Introduction 



The coupling of the vector and pseudoscalar heavy-light mesons to the pion, gH*Hw , 
and the related 7)q defined through 1 

2^/m H m H * _ 

9h*h* = 7 9q, (1) 

Jn 

have been extensively studied in the literature, both for phenomenological and theoretical 
reasons (a list of results can be found in refs. [1-5]). qq is particularly important in the 
studies of the g 2 -dependence of the form factor F + (q 2 ) which encodes the non-perturbative 
QCD dynamics of D — > n and B — > n semileptonic decays. To better illustrate this, we 
consider the dispersion formula for F+^ n (q 2 ) : 



F^ 2 ) = tt~ — + - f dt , (2) 

2ms* 1 <T_ W( m ,+ ms )2 t-q 2 -ie 



f -I , 

9 

mi 



where me* and /b* are the mass and the decay constant of the £?*-meson. Near the 
£>*-pole, F + can be approximated by the first term of the r.h.s of eq. (2), which is why 
9b*Bti is so important. Inversely, from the g 2 -dependence deduced from the lattice data 
for F+^ n , one can extract indirectly the value of Qb*b-k (see e.g. [6]). Moreover the need 
for a more accurate estimate of <tq became recently more evident since the guidance of 
the chiral extrapolations by the predictions of the heavy meson chiral perturbation theory 
necessarily involves g Q [7]. In the effective theory based on the combination of the heavy 
quark and chiral symmetries [8, 9] one of the essential parameters is the coupling ^ which 
appears as the infinite quark mass limit of Qq : 

9 Q = 9oo + O [l/ml) . (3) 

While the B* — > Bn decay is forbidden by the lack of phase space and therefore the 
coupling 9b*b-k cannot be measured in experiment, the coupling 9d*d-k has already been 
measured. CLEO reported recently [10]: 

g D * Dn = 17.9 ±0.3 ±1.9, (4) 
i.e g c cxp - = 0.61 ±0.01 ±0.07, (5) 

where we used eq. (1) to obtain g c from go*D-K- Since the pion emerging from this decay is 
soft, the experimental measurement of T(D* + ) is difficult. The experimental confirmation 
of the CLEO results would be highly welcome. 

Regarding the theoretical estimates the situation is still quite unsatisfactory. Typical 
values for <7q, as obtained by using the QCD sum rules (QCDSR) are low, around 0.3 [1], 
in disagreement with <? c exp '. Such a disagreement with the experimental value was not, 
however, observed with similar couplings, such as g puJn , gNNn, or (fe*A c7r . To overcome 
that problem, in ref. [11] it has been proposed to include the first radial excitations in the 

1 The subscript Q e {c, b} refers to the heavy quark constituent of the heavy-light meson. 
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hadronic part of the sum rules, below the duality threshold. As a result the coupling g c 
becomes much larger (g c ~ 0.7), while the similar modification of the sum rule leaves the 
decay constants (fo and fo*) unchanged. 

Predictions of various quark models are in the range 0.3 < g^ < 0.8 [12]. In particular, 
in the model with the Dirac equation a value of #00 around 0.6 was found before the exper- 
imental result became available [2]. On the lattice, a quenched calculation by UKQCD [13] 
led to (Too = 0.42 (4) (9). This exploratory study was carried out in the static heavy quark 
limit where the difficulties to isolate the ground state are well known. In addition it was 
performed on a coarse lattice, with a low statistics and with only two light quark masses 
from which it was extrapolated to the chiral limit. Finally, the indirect measurements using 
the lattice data for the form factors F + and Fq [6] have extracted rather small values [6]. 

In our recent paper [4], we presented the results of the quenched lattice simulation in 
which we used the propagating heavy quarks with the mass around the physical charm 
quark. We obtained, g c = 0.67 ± O.OSt^Qg, an d observed no apparent dependence on the 
heavy quark mass. Actually the 3 values of Qq measured directly in the neighbourhood of 
the c-quark mass were compatible within our error bars and shew only a very small negative 
slope (see section 4.3 of ref. [4]). This situation allowed a safe interpolation to the c-quark 
but the determination of g^ relied heavily on the hypothesis that the observed behaviour 
was valid over the whole mass-range. In order to check the validity of this asumption we 
have performed a lattice quenched computation of g In the static heavy quark limit, the 
outcome of which is presented in this paper. Our result is 

#00 = 0.48 ±0.03 ±0.11 , (6) 

where the errors are statistical and systematic, respectively. Using this and the results we 
presented in ref. [4], we were able to interpolate to the 6-quark. We get 

g b = 0.58 ±0.06 ±0.10, (7) 
gB*B-n = 47±5±8. 

By comparing eqs. (6) and (7), our results suggest that the I/tub corrections to 'goo are of 
the order of 20 - 30%. 

The paper is organized as follows: after this introduction we explain in Sec. 1 the 
relation between g and the matrix elements that we compute on the lattice; in Sees. 2 
and 3 we present the details of our lattice study and show our main results; Sec. 4 is 
devoted to the discussion of systematic errors; in Sec. 5 we interpolate to the physical 
.B-meson mass; finally we summarize our findings in Sec. 6. 



1 Extracting 

The coupling gwH-w (with H e {D, B}) is defined by the matrix element 

( H{p') n(q) | H*(p, A) ) = g H . Hv (q ■ e\p)) , (8) 

where q = p — p' and e x (p) is the polarization of the vector meson. To determine gH*Hn, 
we can consider the matrix element of the divergence of the light-light axial current, 
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Ay = q^fTf^q. In the limit where the pion is soft, we can write 



(H(p')\q^\H*(p, X))=g H ^ 



q ■ e x (p) 



mi 



x Uml + 



(9) 



where the dots stand for terms suppressed when q 2 is small and close to the m^-pole. To 
bring out Qh^h-k from (9) we express the l.h.s in terms of the form factors, A),i,2 , evaluated 
on the lattice. The standard parametrization reads 



(H(p')\A»\H*(p,\)) = 2m H *A (q 2 f—^-q» + (m H *+m H )A 1 (q 2 ) 



+A 2 (q 2 ) 



e x ■ q 



mi + ran* 



p^+p'^- 



m, 



m 



(10) 



Note that when q — 0, the soft pion limit is verified, since q 2 = (m#* — van) 2 ~ 0. 
Moreover, if we choose q = and p — p ' — 0, eq. (10) simplifies to 



(H\A l \H*) = (m H ,+m H )A 1 (0)e > 



(11) 



As discussed in our previous paper [4], when q 2 = the form factors A),i,2 are related 
to one another because the g 2 -poles in eq. (10) are unphysical and therefore their residues 
have to cancel. This leads to the relation 



2m H *A (0) = (m H * + m^A^O) + (m H * - m H )A 2 (0) . 



(12) 



When taking the divergence of the axial current in (10), only the term in A remains on 
the r.h.s. so that by using eq. (9), we arrive at 



9H*Htt 



2m H *A (0) 



1 



[(m H * + mjr)Ai(O) + (m H * - m H )A 2 (0)} 



(13) 
(14) 



Finally, inserting the definition (1) leads to the expression for gq 
g Q = - — -Ai{0) + — -A 2 {0) 



2^/m H m H * 



(15) 



Owing to the heavy quark symmetry, the vector and the pseudoscalar heavy-light mesons 
are degenerate in the static limit and eq. (15) simply becomes 



?oo=^l(0). 



(16) 



Thus, the value of g^ (in the static heavy quark limit) is simply given by the form factor 
A 1 at q = 0. 
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2 Strategy for the calculation of on the lattice 

To determine the form factor Ai(0) from the lattice evaluation of the matrix element 
(H {A 1 "-] H*), we use eq. (11) where H and H* now refer to the pseudoscalar and vector 
mesons containing an infinitely heavy quark. We consider the two- and three-point correla- 
tion functions, C2 and 63^, with local (L) and smeared (S) interpolating fields for mesons 
consisting of a static heavy and a propagating light quarks. 

In the case in which both the source and the sink are smeared, the two-point functions 
are defined as 



c* s (t x ) = (£ p s (*)p s \o)) = \ E (E Vf(x)Vf\0)) , 



;i7) 



i=i 



where P s and are the smeared interpolating fields of the pseudoscalar and vector 



mesons. Cf L (i,;) is obtained by simply replacing P s — > P 1 



975 



h, V s 



V 1 



where h stands for the static heavy quark field. The equality between the pseudoscalar 
and vector meson correlation functions in eq. (17) is a consequence of the heavy quark 
symmetry. When the excited states are decoupled and the ground state isolated (at large 
enough t x ), the two-point function becomes 



C s 2 s (t x ) = (Z 



where Eh is the binding energy of the static heavy-light meson, and 

z s_ (Oglg) 
(2m H y/ 2 ' 



(18) 



(19) 



The three-point Green functions with smeared interpolating fields are defined in the 
following way, 



0't.*»)=(E P s (y)A l/ (x)vf(0)} 



x,y 



(20) 



o<t x <t v 



where the vector and pseudoscalar mesons are inserted at the origin and at a fixed time 
t y , respectively. The insertion time t x of the (local) axial current (A^ = q^/y^q) varies in 
the range — t y . The matrix element we are considering, (if \A^ \ H*), is then extracted 
from the following ratio: 



R ss (t x ) = 



1 3 



CH(t x )Z s Z s 



3 , =1 C2 S {t x )C2 S {t y — t x ) 

1 * - (H \A t (t x )\ H*} 
2m- 



(21) 



i=l 



= A 1 (0) = g c 



o<.t x <g:ty 



The time dependence cancels in the ratio, so that a plateau in time of R ss (t x ) directly 
leads to the value of Ai(0), which for m q — > becomes the desired coupling ^oo. In what 
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follows we will use the notation in which Ai(0) = f^o, even though this is strictly true only 



for m q — > 0. 



In our simulation we choose to work with the usual (improved) Wilson light quarks 
and with the static heavy quark by using the Eichten-Hill action [15]. In terms of quark 
propagators and the Wilson line, the three-point correlation function (20) can be written 

as 



0**) = <E Tr 



x,y 



1 + 7o o0 



)u, 



(22) 



where S q (x; y; U) = q(x)q(y) is the propagator of the light quark q. In our study the two 
light quarks are degenerate, m u = m d = m q . The brackets (•••)[/ indicate the average over 
the background gauge field configurations ([/), whereas is the ordered product of links 
in the time direction (Wilson line), 



P* = 5(x-y) 11 U t (x + tJ), 



(23) 



tz — tx 



with t being the unit vector in the time direction. 

To improve the statistical quality of the signal we implement the recent proposal of 
ref. [17] and replace the simple link variables in the Wilson line by the so called "fat links", 
defined as 



U t (x) - U?\x) 



i=x,y,z 



Staple 



(x,x + t) + U*f ple (x,x + t) 



(24) 



In order to isolate the ground state from the correlators (22) at the smallest time sepa- 
rations, it is necessary to devise an efficient smearing procedure. We use the proposal of 
ref. [14] and replace the quark fields q(x) by 



q(x) - (r+hMr) 



r=0 



=x,y,z 



Y[Uf(x+{k-l)i) 



q(x+ri) + 



Ft, 



x—ki ) 



q{^x—ri) 



(25) 

The wave function <f>(r) is chosen in such a way that the overlap with the ground state 
is increased. We take 0(r) = e~ r ^ Rb , where Rb is a parameter which is fixed by requiring 
that the smearing be optimal. Note that it is not necessary to normalize the wave function 
since the normalisation factors cancel in the ratio (21). The smearing also includes the so 
called fuzzing: the gauge links U F are fuzzed by N F iterations of the following procedure 



U j (z)^U?(x)=V 



C Uj(x) + J2 [Uf taplc (x, x+j) + U^ix, x + j) 

i+3 



(26) 



where V is used to project the fuzzed fields onto SU(3), and C is a real parameter. 
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3 Conditions and results of the simulation 

The main results presented in this work are obtained from the simulation performed on 
a 24 3 x28 lattice with periodic boundary conditions, at f3 — 6.2. We collect 160 independent 
SU(3) gauge configurations in the quenched approximation. The Wilson fermion action is 
non-perturbatively 0(a) improved with csw — 1-614 [18]. The light quark propagators are 
computed with the following Wilson hopping parameters: 

K q = 0.1344 91 , 0.1348 92 , 0.1351 g3 , (27) 

corresponding to quark masses around the physical strange quark. 

We will present our results using our preferred set of parameters, t y = 13 a, R max = 5 a, 
Rb = 3.0 a, Np = 5 and C — 4. The motivation for this choice of parameters, the effect 
of the smearing (25) and of the modification of the static heavy quark action (24), will be 
discussed in section 4, where we will also quantify the systematic effects involved in our 
calculation. 

3.1 Study of the two-point functions 

Light-light operator: The axial current appearing in eq. (20) consists of two light quarks. 
Even though we use the improved Wilson quarks, for our purpose there is no need to 
improve the bare axial current since the form factor A\ is insensitive to the presence of 
the 0(a) improvement term. On the other hand, the improvement of the renormalization 
constant is of course necessary. We use 

ZM) = zf{gl)(l + ~b A {gl)ap) , (28) 

where the quark mass ap is obtained from the axial Ward identity, d^A^ = 2pP (P = qj^q). 
The non-perturbatively evaluated renormalization constants are: = 0.81(1) [19,20] 

and 6^ = 1.19(6) [19] at /3 — 6.2. We have checked that the inverse lattice spacing obtained 
from the pion decay constant is consistent with the one we reported in our previous paper. 
Since the temporal extension of the lattice used in our previous paper is much larger than 
the present one we will use, whenever needed, a~ 1 (f 7T ) = 2.71(12) GeV from ref. [4]. 

Static heavy-light correlators: The binding energy Eh and the constant Z are obtained 
from the two-point correlation function as indicated in eq. (18). To identify the plateau in 
time where the ground state is properly isolated, we study the effective binding energy, 

The superscript U SS" refers to the situation in which both interpolating fields are smeared. 
Note also that the projection of the static quark onto the positive energy states 
reduces the contamination of the signal coming from the opposite side of the periodic 
lattice. Therefore, the time reversal symmetry over the time dimension of the lattice is not 
preserved in the heavy quark limit and the two-point functions are no longer symmetric 
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Figure 1: Signals for the effective binding energies obtained by using eq. (29) in case of Smeared- 
Smeared ("SS") and Smeared-Local ("SL") sources. In both cases the light quark mass corre- 
sponds to K q = 0.1344, whereas for the static quark we use the "fat" links as discussed in the 
text (see eq. (24)]). 



with respect to t = T/2. We illustrate in fig. 1 the signals for the effective binding 
energies of the Cf and C 3L correlators. By fitting to the form (18) in the time interval in 
which the plateau is observed, we extract Eh and Z. We have also considered a fit with 
two exponentials, but did not observe any noticeable effect. The double exponential fit 
becomes necessary only for when performing simulations without the "fattening" indicated 
in eq. (24) because in that case the plateaus of E eS (t) are much shorter. 

We first calculate the 5-meson decay constant, /^ s atlc . Rather than providing an accu- 
rate value for this quantity (a precision determination of /^ 3 atlc , has already been made in 
ref. [17]), our intention is only to test the validity of our computation. Z L is not obtained 
directly from the local-local correlators because the corresponding correlators are much 
less stable and statistically noisier than the smeared-local ones. For this reason Z L is 
extracted from 

-> Z L . (30) 

In table 1, we present the values of Eff , E H , Z s , and Z L for a static heavy quark and light 
quarks with masses corresponding to the k 9 -'s specified in eq. (27). The decay constant 
ytrtatic j g obtained from Z L using = 5.37 GeV, i.e. the physical B s -meson mass, 



j static = z£ a *Z L J—ar 3/2 , (31) 



Rz(t) 



Z b x 



C^(t) 



a*® 
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aEff aE s H L Z s Z L 


static-gi 

static-g2 
static-g 3 


0.760(4) 0.758(4) 64.8 ± 2.3 0.126(3) 
0.748(5) 0.745(5) 60.8 ± 2.4 0.121(3) 
0.740(5) 0.737(6) 57.8 ± 2.5 0.117(3) 



Table 1: The values of the binding energies obtained by using the Smeared-Smeared and Smeared- 
Local two-point correlation functions. The corresponding constants Z s and Z L are also given (see 
eq. (18)). The time intervals chosen for the fits are: t/a € [8, 13] for Eff , Z s , and t/a G [9, 13] 
forE s H L ,Z L . 

and Z^ at = 0.77(3), the static heavy-light axial current renormalization constant com- 
puted non-perturbatively in ref. [16]. Notice that this value of Zf at is obtained by using 
the Eichten-Hill static heavy quark action, i.e. without including the effects of the fatten- 
ing (24). In ref. [17], however, it has been argued, on the basis of the observed behaviour 
of the step scaling function, that the effect on Z^ tat of the modification of the static quark 
action is only very small. Therefore the value we are using is certainly sufficient to achieve 
our present goal which, we recall, is to check the compatibility of our results with previous 
ones. 

In fig. 2, we provide an illustration of the signal (plateau) of the ratio R z (t), and plot 
the corresponding f^ Uc , as obtained by using eq. (31). From the linear interpolation to 
the strange quark mass (corresponding to m 2 ss = 2m 2 K — m 2 ), we obtain f Bs , which we 
quote as 

y stat 

f*f c = -A- ( 251 ± 20 MeV ). (32) 

This value agrees quite well with the one, f^ tlc = 225 ± 10 MeV, recently obtained also 
at = 6.2, in ref. [17]. 



3.2 Study of the three-point functions 

We observe that, for t x not too close to or t y , the three-point function, Cf 4 f (0, t x , t y ), 
introduced in eq. (20), does not depend on t x . The reason for this is that, in the static 
limit, the ground states of the vector and pseudo-scalar mesons are degenerate. Therefore, 
by writing CH(0,t x ,t 

y) as 

£5>|p(g|# fe > (H k \Mt x )\H;) l^io), (33) 

k I 

where k and / label the excitations, one sees that in the diagonal contribution (k = I) the 
exponential time dependence writes 

^(E^-EH^-E^ty = e -E Hk ty ? ( 34 ) 
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Figure 2: On the left, we illustrate the signal for the ratio R z (t) in the case k = 0.1344. The fit 
of the plateau (t/a G [9; 13] J to a constant gives Z L , which is then used to compute f^ atlc (see 
eq. (31)). The corresponding decay constants are shown in the right plot (empty circles). The 
result of the linear interpolation to the strange light quark, i.e. /^ atlc , is marked by the filled 
circle. Notice that on the right plot only the statistical errors due to Zl are displayed. 



which is independent of t x (t y is kept fixed). The same is true for the ratio (21) in which 
the denominator is composed of two correlators Gi which are degenerate in this case so 
that the exponential t x dependence explicitly cancels out, and one ends up with 

R SS (t x ) = f 2 ii{t J St • (35) 
Z s e~ h H % 

In other words, the plateau for the ratio R ss (t x ) is equivalent to a plateau of Cfjf(tr). 
A remarkable feature is that C3 already develops a plateau at very a small time t x . The 
physical interpretation could be that in the static limit, the non- diagonal matrix elements 
of the axial current simply vanish, 

(F fe |A(4)|^r) MZ ~0, (36) 

which implies that the ^-dependence in the l.h.s. of eq. (34) vanishes as well when k ^ I. 
Note that this is the case in a quark model framework where the axial current reduces to 
a light quark spin operator. Since the spatial wave functions of the ground state and the 
excited states are orthogonal, the spin operator does verify eq. (36) 2 . 

As shown in eq. (21), the value of the bare coupling ~gf£\ is obtained from the fit of 
R ss (t x ) to a constant. We illustrate the signal for R ss (t x ) in fig. 3. The corresponding 
results are listed in table 2, where we also give the final results for = Z I A 'g^ , for our 
set of light quarks (27). We also include the values of when the light quark is either 
the strange or the u/d quark. From the linear extrapolation to the chiral limit 3 , shown 

2 This statement is verified by our results using the static heavy quark. It is further corroborated by the 
matrix elements obtained by using the propagating heavy quark with a mass around the physical charm 
quark: we see that, as the heavy quark mass is increased, the three-point function Cjf^ (0, t x ,t y ) begins to 
develop a plateau for smaller and smaller values of t x . 
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Figure 3: Signal for the ratio R ss (t x ) (see eq. (35)). Illustration is provided for n qi = 0.1344. 
From the fit in t/a G [3, 10], we obtain the value of gj^J . 

in fig. 4, we finally obtain, 



The results presented in this section have been obtained after having fixed t y = 13a < T/2. 
This choice was made after we performed several simulations with t y / a e [8, 13] and studied 
R ss (t x )Zg for each t y . We have checked that the result for R ss (t x )Zg does not change as 
soon as t y gets larger than 10a. To separate the sources as far as possible we have chosen 
to fix t y at 13 a. 

4 Smearing, Fattening and Systematic uncertainties 

As we already mentioned, to improve the statistical quality of the correlation functions 
computed with the static heavy quark action, we used the "fattening" procedure [17]. 
In addition, we use the smearing of the heavy-light interpolating fields [14], which helps 
isolating the lowest bound state at lower time separations. In this section we explain how 
the parameters of the smearing procedure were chosen (25) and show the benefits of using 
the "fat link" static heavy action instead of the standard Eichten-Hill one. After that 
discussion, we enumerate the sources of systematic uncertainties and estimate the related 
errors. 

3 We defer to subsection (4.3) the discussion of the validity of this procedure and of the uncertainties it 
induces. 



0.48(3) . 



(37) 
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-(0) - 


static-gi 

static-g 2 
static-g 3 


0.67(2) 0.58(2) 
0.65(3) 0.55(2) 
0.63(3) 0.53(2) 


static-s 


0.64(3) 0.54(3) 


staXic-u/d 


0.59(3) 0.48(3) 



Table 2: Bare and renormalized values of goo. The time interval chosen to ht for a plateau in 
the ratio R ss (t) is t/a G [3, 10]. 



4.1 Smearing 

The smearing procedure involves four parameters: R max , Rb, Np and C (see eqs. (25- 
26)). The value of R ma x represents half the size of the hypercube where the smearing wave 
function lives and it is set to the maximal value allowed by the memory limitations. In our 
case, R max = 5 a. The parameter Rb, which appears in the wave function 0(r) = e~ r l Rb , is 
not fixed by any physical requirement other than improving the shape of the plateau for 
Rz- The smeared-local correlators are much more sensitive to the variation of Rb than the 
smeared-smeared ones. In fig. 5 we show the shape of Rzif) for Rb G [1.0a , 4.0 a]. We see 
that if Rb > 3.0 a, the value of Rz(t) on the plateau, and thus also the corresponding f^ 10 , 
remains unchanged. The reason is that, as the wave function <f>(r) lives in a typical volume 
of radius R ma x = 5 a, the shape of the normalized wave function term [(r + l/2) 2 0(r)] in 
eq. (25) does not vary much when changing Rb in 3a — 5a. Of course, all the physical 
quantities are built in such a way that, at the end, the effect of the normalization of <j)(r) 
(i.e. the unwanted effect of the smearing) cancels. On the basis of these observations, we 
fix R b = 3a. 4 

The remaining two parameters, C and Np, enter in the fuzzing procedure (26). As far as 
C is concerned, it turns out that its value does not make any significant effect on our result. 
We varied its value in the range C G [0.8 , 4.0], and the shape of Rzif) does not change 
at all (of course, within our present statistical accuracy). On the other hand, the number 
of fuzzing iterations, Np, is chosen to be Np = 5. The effect of each iteration is to fuzz 
the spatial links by the gauge fields in the immediate neighborhood. By several iterations 
we can restore the effect of the wave function inside the volume of radius Rb, where <p{r) 
is supposed to live. As Np iterations cover a volume of typical radius (Np + l) 1//2 a, we 
should require Np ~ 8 to fill a volume of radius -Rb = 3.0 a. It turns out, however, that for 
Np > 6, the signal for Rz begins to increase sharply: a large number of fuzzing iterations 

4 As a side remark, we note that the physical Compton wavelength of the £)-mesons Ad is given by the 
mass splitting between the D and the D\ mesons masses. Its value, \d — 2.2 fm, is somewhat similar to 
Rb = 3.0 a ~ 2 fm, that we use in our simulation. 
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Figure 4: Linear extrapolation ofgoo to the chiral limit. The empty circles correspond to the 
simulated points whereas the rilled symbols are the result of extra(inter)polation. 
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Figure 5: Effect on the ratio Rz(t) of the variation of Rb- For Rj, > 2.0 a the value of Rz remains 
almost unchanged in the region t> 7 a. For this illustration we used the Eichten-Hill static heavy 
quark action, while the light Wilson quark k is 0.1344. 
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destroy the quantum fluctuations as it was observed in the cooling procedures used to 
study instantons [21]. We checked that Rz, and thus /^ tatlc , remain stable for N F 6 [3,5], 
and for the final simulation we chose Np = 5. 



4.2 Effect of "fattening" the static heavy quark propagator 

In this subsection we show that replacing the link variables in the Eichten-Hill action 
by the "fat" ones (see eq. (24)), does make the plateau for the effective binding energies 
much larger, including a visible reduction of the statistical noise. In fact, after replacing 
Ut — > Ul &t in the Wilson line, the plateau gets extended to t ~ 12 a — 13 a, which is actually 
the upper limit (T/2 = 14a). From the comparison of the effective energy, eq. (29), 
obtained by using the "fat" Wilson line and of the standard one (i.e. without fattening), 
we can clearly see the improvement in the signal. This is illustrated in fig. 6, where in 
the case of "fat" Wilson line, the statistics is half what it was in the other case. This 
improvement is still more striking in the case where one uses the so-called Hyp-fattening 
(cf. refs. [22,23]) instead of the standard one (see ref. [24]). We see that not only the 
plateau appears much clearer, but the signal remains longer in time. We also see that the 
two values of the effective binding energies are different, as expected, since the heavy quark 
lattice action was modified. In table 3, we compare the results obtained in both cases, i.e. 
with and without using the fattening procedure. 
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Figure 6: Comparison of the signals for the effective binding energies Egff and Z L when using 
either the usual Eichten-Hill action or its "fattened" version (see eq. (24)). 



Before closing this subsection, we note also that, as far as the stability of the signal is 
concerned, the effect of replacing the standard Eichten-Hill action by the fat-link one is 
much more rewarding in the case of the heavy-light decay constant than for our #00. This 
is so because the operator whose matrix element is related to f^ atlc includes a static quark, 
whereas the one related to ^ does not. The fat link action however helps improving the 
statistical quality of the signal for goo. We note that, besides this improvement, the use 
of fattened links leads to a dramatic decrease of the value of /j* atlc with respect to what 
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Eichten-Hill 'Fattened' 


ystatic [ Me y] 


410(56) 251(20) 
0.58(8) 0.50(3) 



Table 3: Values of f§f c and goo as computed from the usual Eichten-Hill static quark action 
and the new procedure with "fat" links. Parameters of both simulations are kept identical. In 
particular, to determine g^ we have taken in this case t y = 10 a. 

had been found with the Eichten-Hill action. The fact that the plain EH formalism faces 
problems when trying to study fs has been known for some time. In refs. [25, 27] it was 
shown that the smearing led to systematically lower values 5 . We show here that a similar 
effect results from the use of fattened links. In ref. [24] we show that the Hyp-fattening 
confirms fully the present results. 

4.3 Systematic errors 

• Smearing: 

We combine in quadrature the effects of the variation of all the parameters of the 
smearing procedure on the resulting value for (see discussion in sec. 4.1). We 
estimate that systematics to be ±15%. 

• Discretization errors: 

In our study we implemented the full O(a) improvement of the Wilson QCD action 
and the axial current. As we discussed in the text, the improvement of the bare 
axial current does not influence the value of g^. As for the renormalization constant, 
we used the non-perturbatively determined value, including the coefficient 6^, which 
ensures the elimination of the artifacts of 0(ap). 

Our main result is obtained from the simulation at f3 = 6.2. In order to study the 
0(a) effects, we also performed the simulations at f3 — 6.1 and (3 = 6.0, keeping the 
physical volume approximately the same and rescaling the smearing parameters. In 
fig. 7, we show the chiral behavior of goo, as computed in all three simulations. We 
see that both the results accessed directly from our simulations and those obtained 
after the linear extrapolation agree within the error bars. To give the reader a more 
quantitative insight into this effect, we collect in table 4 the results for 'g OQ and /j^ atlc 
from all three simulations. 

Since in this work we do not deal with propagating heavy quarks (which are usually 
the dominant source of discretization errors in the heavy-light meson observables) , 
and since we do not observe any noticeable discretization error on from our 

5 We thank the anonymous referee for drawing those references to our attention. 
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Figure 7: Comparison of the chiral extrapolations of from three simulations at (3 = 6.0, 
6.1 and 6.2. The values extracted from these simulations are compatible within one standard 
deviation. 



simulations at three different lattice spacings, we will not account for any extra 
discretization error in our overall systematic uncertainty 6 . 

• Finite volume: 

In this work we did not carry out a detailed study of the finite volume effects. Instead, 
we will rely on our previous study with the propagating heavy quark where the finite 
volume effects are estimated to be ±6% [4]. This uncertainty will be added in our 
overall systematic error estimate. 

• Chiral extrapolations: 

Since we have only three values for the bare light-quark masses, we cannot make a 
detailed study of the effect of the chiral extrapolation like the one we did in [4] . The 
values of ^ that are directly accessed from our simulations follow a very smooth 
linear behavior when the value of the light quark mass is changed. The result (37) 
is obtained from the linear extrapolation. To that we will add an error of ±15%, 
as estimated in our previous paper [4]. This error is considered symmetric for the 
reasons explained in detail in ref . [28] : the chiral logs in the full (unquenched) theory 
have a tendency to lower the value of the coupling g^, whereas those in the quenched 
theory do the opposite (increase the value goo). 

6 Notice also that goo is dimensionless so that no discretization errors due to the conversion from the 
lattice to physical units arise. 
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6.0 6.1 6.2 


f§f c [MeV] 

9oo 


250(35) 257(27) 251(20) 
0.45(4) 0.51(3) 0.48(3) 



Table 4: Values of /^ s atlc and 'g 00 from three simulations at (3 = 6.0, 6.1 and 6.2. Consistent 
results are found when comparing these studies. 



After combining the above sources of systematic uncertainty in quadrature, we end up 
with 22% of error to the value given in eq. (37), i.e. 

^oo = 0.48 ±0.03 ±0.11 . (38) 



5 Heavy mass interpolation to the 6-quark 

We can now combine our static heavy quark result, f^o, with those obtained in our 
previous study in which propagating heavy quarks with masses around the physical charm 
quark were used [4], to interpolate to the 6-quark sector. For that purpose we use the 
spin-averaged mass of the heavy-light mesons, i.e. 

_ 3m v + m P 

m H = . (39) 

Motivated by the heavy quark symmetry, we fit our results to the linear and quadratic 
forms: 

9Q = 9oc + — , 40 

^ ^ bi 6 2 / A , x 

9q = 9oo + — + =2- • (41) 

These fits are shown in fig. 8, from which we read off the values at l/m B = (5.314 GeV) -1 , 
i.e. g b : 

g b ^ = 0.55(4)(9), g b (quad) = 0.60(7) (9) . (42) 

As our final estimate we will take the average of the two, and add the difference in the 
overall systematics. Thus we have 

g b = 0.58 ±0.06 ±0.10, (43) 

which by means of eq. (1) gives 

<? B * B7r = 47±5±8. (44) 
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around the charm (empty circles) and our new result at the static limit (hlled circle). The linear 
and quadratic fits defined in eqs. (40) and (41) give compatible results for % (filled square and 
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6 Summary and conclusions 

In this paper we made a study of goo, the coupling of the lowest-lying doublet of heavy- 
light mesons to a charged pion, in the infinitely heavy quark limit. The static heavy quark 
action used in this work is the Eichten and Hill one with fat-links, which ensures a better 
statistical quality of the correlation functions computed on the lattice. Our result, in the 
quenched approximation, is 

<?oo = 0.48 ±0.03 ±0.11, (45) 

where the errors are statistical and systematic respectively. The most significant sources 
of systematic uncertainty are those related to the chiral extrapolation and to the smearing 
procedures used to suppress the contribution of the higher excited states in the relevant 
correlation functions. 

That result is then combined with the ones reported in our previous work [4] , in which 
we used propagating heavy quarks with masses close to the one of the physical charm 
quark. Interpolation to the 6-quark sector leads to 

g b = 0.58 ±0.06 ±0.10 and # B > Bw = 47 ± 5 ± 8 . (46) 

The quantities we have studied in this paper, and qb* Bn, are very interesting because 

1 they are of high phenomenological relevance, 
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2 they escape direct measurement, 

3 their theoretical determinations are very widespread. 

Still, as the attentive reader will have noticed, the errors remain larger large. Referring for 
instance to eq. (28) above, the statistical error is 6% and the systematic one 22%. 

Concerning the statistical error a first factor is obviously the performance level of the 
available computing facilities. Since this is getting higher and higher rather fast and, at 
the same time, new computational techniques have been discovered which greatly improve 
the signal/noise ratio this source of uncertainty will probably be tamed very soon. Some 
of the systematic errors are no serious problems : most of the various renormalization 
and improvement constants we need have now been non perturbatively determined with 
an accuracy at the %-level. Note also that some of the (less well-known) heavy-light 
constants (e.g. Z^ Uc ) simply drop out of our computation. Some other ones are more 
serious : 

The chiral extrapolation (which we have very roughly estimated to contribute by 
15%) can be attacked on both the theoretical side (determination of the chiral logs) 
and the numerical one. 

In this study we have used the static limit only as a limit point for our extrapolations 
from the finite-mass region. It is possible to go beyond this approach and to determine 
higher order -rj- coefficients, which would allow a kind of "educated fit" and to 
decresase the theoretical uncertainties. This is part of an ambitious program which 
has been undertaken by the Alpha-group 3 years ago and is still under progress [29- 
32]. 

There remains the question of the unquenching. Its possibility while using the fatten- 
ing procedures has been established (see ref. [33]) at least in the staggered- fermions 
case and probably with Wilson fermions too. This is one of the tasks we should 
undertake in the future. 

These are just a few hints regarding what will have to be done now but it is impossible at 
present to outline any really sensible quantitative perspective regarding the time evolution 
in the near future of the various uncertainty factors we have just mentionned. Still, pursuing 
in the direction of a "pure QCD" determination of ~g is clearly of primary importance. 
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